Level statistics inside the core of a superconductive vortex 
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Microscopic theory of the type of Efetov's supermatrix sigma-model is constructed for the low- 
■ lying electron states in a mixed superconductive-normal system with disorder. The developed tech- 

nique is used for the study of the localized states in the core of a vortex in a moderately clean 
superconductor (1/A <C t ^ lo^^ = Ep/A'^). At sufficiently low energies e ^ LOrh, the energy level 
statistics is described by the "zero-dimensional" limit of this supermatrix theory, with the effective 
"Thouless energy" uiTh ~ (loo/t)^^^. Within this energy range the result for the density of states 
is equivalent to that obtained within Altland-Zirnbauer random matrix model of class C. Nonzero 
modes of the sigma-model increase the mean interlevel distance u>o by the relative amount of the 
order of [21n(l/LJor)]-^ 



. There is a great deal of activity during the last 1.5 decades directed to the study of electron energy levels and 
D ' wavefunctions in disordered normal metals 0], where they govern low-temperature transport properties. In (s-wave) 
^ , superconductors, disorder is usually of less importance, since the excitation spectrum is gapful and single-electron 
states are almost empty at T ^ A. The situation is quite different in mixed superconductive-normal systems (for 
d . recent reviews see |Q) where the gap in the excitation spectrum can be: i) very low compared to the bulk A, or ii) just 
zero. The example of the first case is presented by an S-N-S sandwich with the thickness of the N region much 
longer than the superconductive coherence length ^ and the mean free path I. Then at generic values of the phase 
difference ip between superconductors the gap in the electron spectrum in the N region is of the order of the Thouless 



energy Eth = D/Lj^ <C A. To calculate the density of states (DoS) fluctuations at e > Eth, and other mesoscopic 
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Q ■ eiieigy ^Th = ^1 

Q ' effects in such systems, Altland, Simons and Taras-Semchuk developed recently |^ a field theory which is a version 
of Efetov's supersymmetric matrix sigma-model explicitly taking into account the superconductive coherence induced 
' in the N metal due to the proximity effect. At low enough energies the effective theory they derived is equivalent to 
^ . the Orthogonal random matrix ensemble at (/? = 0, and to the Unitary ensemble at 3> 5/ Exh (where 5 is the level 
\G ' spacing in the N) . Like in the usual normal- metal case , mesoscopic fluctuations reveal themselves in the quantities 
containing products of both retarded and advanced Green functions, whereas single, say, retarded Green function is 
' determined by the standard equations of the quasiclassical theory of superconductivity. 

, Qualitatively different theoretical problem comes about in the second case, ii), mentioned above, which is realized, 

■ e.g. in the same S-N-S sandwich at (/? = tt or in a variety of situations where the phase of the order parameter rotates 
0^ ■ due to the presence of an external magnetic field. Now the DoS is non-zero at arbitrary low energies, and quantum 

interference due to Andreev scattering strongly affects even the average DoS (p(e)) at e ~ (5. General approach to 
this kind of systems, characterized by the zero averaged over the whole system value of the superconductive order 
parameter, was initiated by Altland and Zirnbauer (AZ) j5|, who employed a generalized random-matrix approach. 
They have shown that the particle-hole symmetry of the Bogolyubov-De Gennes (BdG) Hamiltonian leads to important 
^ ■ constraints to be imposed on the random-matrix Hamiltonians. Precise form of the constraint depends on the presence 
O ■ or absence of: a) time inversion symmetry, and b) spin rotation symmetry. Thus AZ identified 4 additional (with 
^ ' respect to the standard Wigner-Dyson theory) classes of random-matrix ensembles appropriate for the description of 
J> . mesoscopic fluctuations in this kind of S-N-S systems. Crossover between such classes has been considered in using 
the space-independent supermatrix sigma-model. While the AZ approach is certainly highly suggestive, it has the 
same general limitation as any ad hoc random matrix theory, i.e. the limits of its applicability to some real physical 

■ system are left undetermined. Therefore we consider it highly desirable to develop a fully microscopic field theory for 
the mesoscopic fluctuations in S-N-S systems without a spectral gap. 

In the present Letter we develop a microscopic field-theory approach to the generic example of the ii) type of 
systems, namely, to the core of a superconductive vortex. It is known since Caroli, De Gennes and Matricon (CdGM) 
that the BdG equations for the electron states near the Abrikosov vortex possess localized solutions with energies 
well below the bulk A. The spacing between these localized levels, wq, is of the order of A/ {kpS,) and disappears in the 
quasiclassical limit kp^ oo. Thus it was tempting to consider the vortex core as a kind of a "normal tube" inside 
a superconductor and, indeed, in many cases such a simplified picture was found Q to be at least qualitatively 
correct. Later it was demonstrated jlG] that the presence of a quasi-continuum spectrum branch localized on the 
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vortex follows from general topological arguments; actually, the number of such "chiral" branches coincides with the 
topological charge of the vortex. However, it is not always possible to consider the chiral branch as a continuos one, 
as it stands in the quasiclassical calculations. It was shown recently |l^Jl^, that the discreteness of the localized 
energy levels becomes of real importance in layered superconductors at sufficiently low temperatures. In the previous 
paper |ll| we employed the AZ phenomenological approach to find low-current nonlinearities in the current-voltage 
relation in a mixed state of a moderately clean superconductor (the mean free path I ^, but / ^ ^{kpO)- such 
a case the inverse elastic scattering time l/r is much larger than inter level spacing loq, therefore the applicability of 
an appropriate random-matrix model (which is, in fact, class C of the AZ classification) seems quite natural. Vortex- 
induced dissipation in another limiting case of a super-clean superconductor with extremely low concentration of 
impurities (/ ^ kpi"^) was considered recently by Larkin and co-workers [Q. Here the system of electron levels in the 
core was found to be extremely strongly correlated and almost integrable, with the properties qualitatively different 
from the moderately clean case In the present Letter we again consider a moderately clean limit wq ^ 1/t <C A, 
now within microscopic approach starting from the BdG equations in the presence of impurities-induced Gaussian 
random potential. We derive the conditions under which the AZ class C statistics is indeed realized in the vortex core, 
and estimate the scale of non- universal corrections to it. Throughout all the paper we consider purely 2-dimensional 
superconductor, which is a good approximation for the case of sufficiently strong layered anisotropy; more quantitative 
conditions can be found in |]lT| . 

Below we present a brief description of our method and results, whereas their detailed presentation is postponed 
to a forecoming paper . It was already mentioned above that in the present problem even the calculation of the 
average single-particle quantities is not trivial and cannot be done within the quasiclassical theory, as long as low 
energies comparable to the level spacing are considered. Thus our goal here is to derive a field-theory technique 
for the calculation of the average DoS (p(e)) = S{e — ej)). To average the Green function over disorder, we use a 
standard trick of representing it as the functional integral over both Grassmann (x) and usual complex (S) fields 
which combine into the superfield The most direct way would be to work with real-space-dependent superfield 
^(r) corresponding to the standard representation of the BdG Hamiltonian. In this way we would obtain a field 
theory in terms of (3(r) supermatrix, depending on two spatial coordinates rx^Vy. On the other hand, low- lying states 
of the chiral branch depend upon a single quantum number only (in the absence of disorder it is just the angular 
momentum), as well as for a generic ID problem. Therefore, in the basis of such states the BdG Hamiltonian can be 
represented as a random N x N Hermitean matrix (where N ~ A/wq is the total number of the localized states in the 
core) of the certain structure and symmetry which we will discuss below. In the clean limit, l/r ^ A, the admixture 
of delocalized (e > A) states to the low-lying ones can be neglected. Thus it is convenient first to reduce the full 2D 
problem to a sort of random matrix problem that can be further reduced to the ID field theory, explicitly containing 
the chiral spectrum branch only. 

In the basis of the chiral CdGM states *^(r) = A( J^,_i/2(fcFr), J^+i/2(fcFr))'^e^''''e"^('') determined in (here 
A ~ ^Jkp/S, is the normalization constant, is the azimuthal angle in the real space, /x e [— A^/2, N/2] is the angular 
momentum that takes half-integer values, and K{r) = (l/hvp) /o*^ A(r')(ir'), the full Hamiltonian takes the form 
{l2\H\fi') = ujo^Sfj_,fj_' + where the second term is due to the random white-noise impurity potential U(r) with 

the variance {U{r)U{r')) — 5{r — r')/ {2tti/t); correspondingly, in the functional integral one should use /i-dependent 
supervector <I>^ instead of the superfield $(r). This Hamiltonian obeys the symmetry 

H = -jH^j^; (^|7|A*') = (-ir+'/'W> (1) 

which follows from an identity 5'_^(r) — (— l)''+iiTj,^* (r) that reflects the basic symmetry property of the BdG 
Hamiltonian. Here we introduced the Pauli matrices T^^y^z in the 2x2 Nambu space. 

The standard way to solve a complicated random matrix problem is to represent it in a form of the effective field 
theory. In order to reduce the random matrix problem given by Eq. (^ to the ID field theory we make a continuous 
Fourier transform (considering N as very large, which is possible since the energy range e <C A is studied below) from 
the momentum variable fi to the "angle" (j) G [0, 27r), so our superfield will be defined as <&((/') — 'J'a'^"*^^"^^^'''^- 
Now we can write down an expression for the 'partition function' (e+ = e + iS): 

Z^ie) =I^^P^[I (^+ ~ *^o^ - y ) m ^$*(0)F(<^, 0')'i>(0')} D<^*{cl^)Dm- (2) 

Matrix elements V{(j), <j)') of the random potential in the ^i-space obey the symmetry relationship that follows from 
Eq. (|l|) and are given by 

1/(0,0') = -e»('^-'^')-l/*(0 + 7r,,^' + 7r) =yl2 j £vw^^,{r,6)U{T)e-^'^^'-\ (3) 
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The function w^^' can be presented, using summation formulae for the Bessel functions, as 

w^^>{r,0) = (1 - e'('^~'^')) exp |-i/cF''[(sin0 - sine/)') cos 6* + (cos(/) - cos(/)') sin6'] |. (4) 

All features of the theory are encoded in the pair correlator Wicpi, (j)2, (ps, 4>4) — {V {(pi , (j)2)V {(j)^ , (j)4)) where the 
averaging is performed over the Gaussian distribution of the impurity's potential C/(r). Since the typical value of 
kpr ^ kp^, ^ 1, the correlator yV(0i, 4>2, ^3, ^4) is essentially non-zero only when the oscillating exponents in Eq. (Q) 
nearly cancel each other, i.e. when its arguments (pi are nearly pair- wise coinciding p3[ |: 

2 

(/.2, (/-a, M = — T(0i - (/.2 + 7r)(27r)2 U(</,i - 04)^(02 - </)3) - e*(^^-^^)5(0i - <^3 + ^)'5(^2 - </'4 + 7r)l , (5) 

TT L J 

where g = 2A^/'iTVTuo'^k^p ^ l/LUt^T 3> 1 and the kernel T is given by 



cot 



if|'^l>:^ 



iv, if 101 < 
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The (5-function approximation (|^) for the correlator W is valid as long as the scale of the angular variations of the field 
$(</i) (below it will be seen to be £ = [g hi{N / g)\~^) is longer than the actual ||T^ width w{(f>i~(f>2) ~ \N sin{(j)i — (j)2)\^^ 
of those (5- functions. Thus the following derivation is strictly valid under the condition w{£) <ti £ which is equivalent 
to: 



Ty/uJoA > In At. (7) 

Note that this condition is stronger than just the clean limit condition rA 1. In what follows we will assume this 
condition to be always fulfilled. 

The next step of the cr- model derivation is to average the partition function (||) using Eqs. (||), (^. Before doing 
that we need to take explicitly into account the symmetry (|l|), (|^) which amount to the doubling of the number of 
components of the supervector (a similar procedure in the standard approach is related with the time-reversal 
symmetry). Thus we introduce (cf. with a similar procedure in |^) an additional 2x2 "particle-hole" (PH) space 
and define a 4-dimensional supervector field ip{(l)) = 2~^/^($((/)), e*"^<i>*(0 -I- tt))"^. Next we define the bar-conjugated 
supcrfield as = ■4'~^(Tz = [C{(l))tp{(j) + n)]'^ with C(^) = — e~"^(T^Co, where CTz is the Pauli matrix acting in the 
PH space, and 

with FB meaning the Fermi-Bose space. After the averaging over disorder the effective action -4{'!/'} entering the 
functional integral for the retarded Green function ^''^(e) = — « / <i>a*&a exp[^{^/;}]2?'0*I?-0 (where a implies either 
bosonic or fermionic component) can be written as 



Am = ^ J ^e+a. - ^c.o^ - y j m ^ J J ^i^^^'^i " + ^)V;(0i)V'(</'2)V;(</'2)V'(</'i)- (9) 

The second term in the action is very similar to that of the ID tight-binding model with off-diagonal random 
matrix elements with variance decaying as l/|a;|, as long as we are interested in the scale |a;| = \(f)i — 02 -f 7r| <C tt. 
Therefore, we expect the usual ID localization to be absent in our problem just because of the long-range nature of 
the off-diagonal disorder (cf. [|^). 

There is also another suggestive way of considering this term, which helps to gain some intuition about its 
effect. Namely, one can think of the variable (p as an angle associated with the 2D quasiparticle momentum 
p = kp {cos (p,aiii(p}. Then the second term in Eq. (^ corresponds to a 2D particle-hole scattering which is 
strongly enchanced in the forward direction. For such a singular scattering one has to define two scattering 
lengths £ and £tr S> £ (cf. with a similar situation discussed in ||l5|): l/£ (x g J d(p (y{(p) — ghi{N/g), and 
l/£tr o<: g J d(p cr{(p) (1 — cos0) = (7^1, where (j{(p) is the differential cross-section and (p = (pi — (p2 + ■ For 
the careful evaluation of the logarithmically divergent scattering rate l/£, one should use the self-consistent Born 
approximation (SCBA) which takes into account both terms in Eq. (j^). It is equivalent to taking into account of 
the "non-crossing" diagrams that can be generated by a perturbative expansion of exp[yl{?/'}] in powers of g. For 
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e/wo <C l/£, the "crossing diagrams" of the same order in g turn out to be small by the parameter = 1/ ln(iV/g). 
It stands for the usual quasiclassical parameter {kp£tr)^~''' in this effectively ID problem. 

The existence of the small papameter l/\n.{N/g) — 1/ In Ar that allows to neglect the "crossing diagrams" , implies 
that one can derive an effective field theory (nonlinear sigma-model) which describes the low-energy behavior of the 
averaged Green function G^(e) for e/ujo <^ l/£ = g\ii{N/g). This can be done in a standard way by the Hubbard- 
Stratonovich decoupling of the quartic term in Eq. and a further saddle point approximation controlled by the 
parameter \/\Yi{N/g). Because of the symmetry relationship (|l]), (^) and the corresponding relationship between 
and ■0) one has to perform both the local decoupling containing P(0) "0(0) ® '4'{4') a-^id the non-local one containing 
R{(j)i,(j)2) il){(t>i)®i^'^{(t'2)- Under the condition I :$> 1/VN given by Eq. (|^), both decouplings are important in order 
to obtain a correct form for the imaginary part of the Green function in a saddle-point approximation P(0) = Pq, 
R{(j)i, 4>2) = -Ro(0i — 02) which is equivalent to the SCBA: 

Ge(0) = -— a.0(-a,0)e-l*l/V^^^ Po = ^a., i?o(0) - -T(0)G,(-</)), (10) 
where — g\n{N/g), and Tq = Jq^ T{<j))^ « ilniV. In general, m-th Fourier harmonics of the kernel r(0) is 



given by T„ « In ]^ for K rn < VN. 

Mesoscopic fluctuations are known [n| to be described by the slow rotations of the saddle-point solution, which are 
represented in our case as P((/)) = U {(j)) PqU {(f)) , R{(f),4)') — U~^{(j))Ro{4> — (1>')U{(I)'). The corresponding action 
that describes the low-energy spectral properties of the CdGM levels, reads: 

X[Q,C/] = -f T„^//|^r-(0r-02+.)StrQ(0OQ(02)-f / ^ Str (-la.g(0) - ^..L/(0)^^ 

(11) 

where Q{4') = U~^{(j))azU{(f)), and U{(f)) is a 7r-periodic, pseudo-unitary {U~^{(p) = U{(j))) matrix. The action (|ll| ) is 
valid for the energies e <^ / £ — In Ar. 

The supermatrix Q can be represented in the form (5(0) — <7z[l + 1^(0) + \W^{(j)) + 0{W^)\ with the supermatrix 
W being purely off-diagonal in the PH space. Then the symmetry Q = Q and convergence arguments lead to the 
following form for the Wph and Whp blocks: 

Here z is a complex number and Ui are Grassmann numbers. Expanding over VF(0), we obtain in the quadratic 
approximation 

A2[WJ, = J Str^ |2g 1^ |] j + i [ma, - ^) | W^2mM^-2™, (13) 

where Wm is the m-th harmonics of the field 14^(0). Note that in Eq. ( p^ only even harmonics enter; odd harmonics, 
as well as the "longitudinal" modes, have a larger gap of the order oi u)q/£ and are excluded from the sigma-model 
action. 

Eq. ([l3|) sets a characteristic scale L for the angular variations of matrices U (0) : 

l/L = g\ng. (14) 

This scale should be larger than the scattering length £. Only in this case one can restrict oneself by the lowest term of 
the gradient expansion in powers of dU/dcj) that has been used in deriving Eq. (pT]). Comparing to £ = [g\n{N / g)]^^ 
we see that the parameter of the gradient expansion, £/L — In ln(iV/5), is small if the condition (0) is fulfilled. 

The length L determines the angular size of the elementary propagator corresponding to the sigma-model ([T]]). 
In this respect it is analogous to the phase-breaking length for Cooperons (or the system size) in the usual weak- 
localization problem. The fact that £/L <^ 1 in our problem tells us that the problem is essentially not ballistic, 
though it is not diffusive either, since £tr /L = \ng > 1. 

An important property of the action ( pT| ) is that it takes a universal form if U is independent of 0. It is clear from 
Eq. (^|) that at low energies the main contribution to the functional integral comes from the zero harmonics of Q(0), 
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i.e. the problem reduces to the zero-dimensional supermatrix a-model. The uniform supermatrix Q is parametrized 
by 2 real variables (one of which appears to be cyclic) and 2 Grassmann variables, so the final expression for the 
average DoS is 

{P{e))^^^ r de I d,dC-^ [d + cos 6) + 2,C(1 - COS 9)] e-^^^—o^ ^ ^ f 1 - ^^f^fo) \ 
Aluo Jo J l-cosfc''- luq \ 27re/cjo / 

The functional form of this result coincides with the result of the AZ phcnomenological approach However, it is 
expressed via the renormalized mean level spacing ujq ^ ujq (^1 + 2ing ^ • renormalization is due to the contribution 
of higher Wm>2 modes which lead to the decrease of the DoS in the energy range e < luq/L = gujQlng by the 
relative amount of 6luo/ujq = l/(21n(7) <C 1. At higher energies, the correction decreases as Sluq/ujq oc {guoliig)'^ /e^. 
This correction can be found using a general approach , in which the pcrturbative treatment of the non-zero modes 
leads to the "induced" terms in the OD action. It is given by the single- cooperon diagram which is absent in usual 
normal- metal problems | l|,[l6|] ; from the formal point of view, the difference stems from the absence of the BB block 



in the parametrization (12|). A usual p6[ two-cooperon diagram leads to the "induced" term oc {e/uJo)'^{glng)' 
the effective action of the above OD ct- model (cf. ^ where similar result is mentioned) . The possibility to neglect this 
term determines the upper limit of energies where purely OD description is valid: 



e < WT/i = ujoy/g\ng. (16) 

To conclude, we have derived microscopically the supersymmetric field theory for the statistics of the localized 
electron levels inside the vortex in a moderately clean superconductor. Our supermatrix tr-model, Eq. ( pi] ) was 
derived in the main order in the quasiclassical parameter 1/lnAr. Previously proposed random-matrix approach 
1^ is shown to be valid in the low energy range e < ajxh — [('j-'o/''') lii(l/'^o''')]^^^ where zero-dimensional cr-model 
is applicable. Mixing between zero- and higher modes leads to the decrease of the DoS by the relative amount of 
[2 ln(l/cjoT")]~"^ at the energies e < Iii{1/ujot). 
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